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Abstract 

Plasma outflows from gamma-ray bursts (GRB), pulsar winds, relativistic jets, and ultra-intense 
laser targets radiate high energy photons. However, radiation damping is ignored in conventional 
PIC simulations. In this letter, we study the radiation damping effect on particle acceleration via 
Poynting fluxes in two-and-half-dimensional particle-in-cell (PIC) plasma simulation of electron- 
positron plasmas. Radiation damping force is self-consistently calculated for each particle and 
reduces the net acceleration force. The emitted radiation is peaked within a few degrees from the 
direction of Poynting flux and strongly linear-polarized. 

PACS numbers: PACS 52.65.-y, 52.65.Rr, 52.30.-q 
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When charged particles suffer extreme acceleration, radiation loss and damping become 
important in the plasma energetics and dynamics. This is especially true in the ultra- 
relativistic regime, where radiation damping can severely limit individual particle acceler- 
ation. However, conventional Particle-in-Cell (PIC) simulations of collisionless 
plasmas have not included radiation effects. 

One of such an example is high-energy astrophysical phenomena, such as pulsars, blazars 
and gamma-ray burst (GRB) emissions. There are two competing paradigms for the origin 
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of the prompt GRB emissions: hydrodynamic internal shocksp, la III M li| versus Poynting 
fluxes | lflj]. Both pictures require the rapid and efficient acceleration of nonthermal electrons 
to high Lorentz factors in moderate magnetic fields to radiate gamma-rays. In the Poynting 
flux scenario, long-wavelength electromagnetic (EM) energy can be directly converted into 
gamma-rays using the electrons or electron-positron pairs as radiating agents. Such Poynting 
flux may originate as hoop-stress-supported magnetic jets driven by strongly magnetized 
accretion onto a nascent blackhole, or as transient millisecond magnetar wind, in a collapsar 
event 

111] 

or in the merger of two strongly magnetized compact objectsjl^. 
The recent discovery of the diamagnetic relativistic pulse accelerator (DRPA) 

in which intense EM pulses imbedded inside an overdense plasma (EM wavelength 



A ^> plasma skin depth cj u pe ) capture and accelerate surface particles via sustained in-phase 
Lorentz forces when the EM pulses try to escape from the plasma, is particularly relevant 
to the Poynting flux scenario of GRBs. Liang & Nishimura |l4j recently showed that DRPA 
reproduces from first-principles many of the unique features of GRB pulse profiles, spectra 
and spectral evolution. DRPA-like Poynting flux acceleration may occur when a Poynting 
jet head emerges from the surface of a collapsar, or when a new born magnetar wind blows 
out the progenitor envelope. 

In this article we report PIC simulation results of particle acceleration driven by EM- 
dominated outflow (Poynting flux), using a newly developed 2-1/2-D code that includes 
self-consistent radiation damping. We compute the radiation-damped plasma and field evo- 
lutions, as well as observable radiation output, which may be applicable to both laser ex- 
periments and astronomical phenomena. For example, focusing ultra-intense short pulse 
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lasers 1 181] into small focal spots, the intensity may exceed 10 22 W/cm 2 
Schwinger limit, 10 29 W/cm 2 with couner-propagating pulsesj^. Under such conditions, 
electrons can be accelerated to relativistic speed, and the effect of the radiation damping 
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decreases the efficiency of the laser heating 22|. In future, we will apply the code to such a 
case. 

High frequency radiation with A C Ad = c/u pe cannot be included as part of the self- 
consistent EM fields in PIC simulations since the Maxwell-solver cannot capture gradients 
> A^ 1 = uj pe /c, where uo pe = Aime 2 /m e is the electron plasma frequency. 

Relativistic particles can emit up to the critical frequency u c = 37 2 f2 ce , where 7 = 
E/m e c 2 = — v 2 /c 2 is the Lorentz factor and f2 ce = eB/(m e c) is the electron gyro- 

frequency. The ratio of the critical radiation wavelength A c to Ad is given by A c /Ad = 
(2iiu pe ) I \3j 3 Q ce ) , which is 1 because u pe /Q ce < 0.1 in EM-dominated cases and 7> 1. 

To account for high-frequency radiation, we introduce a radiation damping force in the 



24, 



form of the Dirac-Lorentz equation 
given by (see Ref. 4 for detailed calculations 
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251 ] . The classical damping force term f rad is 
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where v is the velocity, and E and B are the self-consistent electric and magnetic fields. 
Here we introduce a non-dimensional factor k rad given by 
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where r e = e 2 /(mc 2 ) is the classical electron radius. 

The first square bracket term of Eq. (0) represents the radiation damping due to the 
ponderomotive force acceleration. The third square bracket term is Compton scattering 
by large scale (A> Ad) electromagnetic field which reduces to Thomson scattering in the 
classical limit 26j|. We note here that Compton scattering with high frequency radiation 
(A Ad) is not included in Eq. (JTJ). 

We restrict ourselves to stay below the quantum- limit, Ml ce < m e c 2 or B < 4.4 x 10 13 
gauss, which corresponds to k rad =< 7.2 x 10~ 3 , otherwise formula (0) fails. This is the 
case for magnetars {k rad ~ 1CT 2 ). We choose k ra d from zero to 10~ 3 in the simulation 



to maximize the radiation effect and run our simulations until r S i m > lO 4 ^^ 1 so that 

\frad\ T sim — \E ex t\£l ce ■ — 

We use the explicit leap-frog method for time advancing Spacial grids for the fields 
are uniform in both x and z directions, Ax = Az — Xd- The simulation domain in the x—z 
plane is —L x /2 < x < L x /2 and < z < L z with a doubly periodic boundary condition in 
both directions. 

As an example to highlight the effect of radiation by intense EM pulse, we here repeat 
the non-radiative simulations by Liang et. al. using the similar initial condition in 
the following simulations. The initial plasma is uniformly distributed at the center of the 
simulation box, — 6 Ax < x < 6 Ax and < z < L z . The background uniform magnetic field 
B = (0, B , 0) is applied only inside the plasma, so that the magnetic field freely expands 
toward the vacuum regions. We choose L x to be long enough so that plasma and EM wave 
never hit the boundaries in the x direction within the simulation time. The initial plasma 
is assumed to be a spatially uniform relativistic Maxwellian, ksT e = ksT p = lMeV, where 
the subscripts e and p refer to electrons and positrons. 

The radiation damping force Q is calculated self-consistently and fully-explicitly as fol- 
lows. The velocity of each particle is updated each single time step At from v(t — At/2) to 
v(t + At/2), using the electromagnetic field and f rad at time t. All the terms in the equation 
(JTJ) can be calculated by the ordinary leap-frog field solver [l|, 0, [J except the displacement 
current dE/dt, which is not given at time t. To calculate dE/dt, particle position x need to 
be updated from t — At/2 to t using velocity v(t — At/2). Next, temporal current Jtif) is 
calculated from v(t — At/2) and x(t). Jt has the first order accuracy in time. Finally, the 
displacement current is calculated using the Maxwell equation, and the radiation damping 
force is calculated for each particle via Eq. (JTJ). To update the velocity, we sum f rad and 
the external force n e t ^oele rati o n . t „da PP l y the Boris rotationQ. 

We run simulations for six different sets of parameters shown in Table HJ k ra d = 0, 10 -4 , 
10~ 3 and Up e /Q ce = 0.1, 0.01. We call Run A-C the weak magnetic field case and Run D-F 
the strong magnetic field case. 

First, we check total energy conservation for Run D,E and F in Fig. In the radiative 
(RD) cases, the energy loss by the radiation E ra d is obtained from the time integral 

E rad = [ \J2 V ' frad) dt. (3) 



In the RD case, sum of E^ in and Ef ie does not conserve (Line 4), but sum of Eun-, Efi e and 
the radiation energy E ra d (Line 3) conserves (Line 5), indicating that the radiation damping 
force is self-consistently calculated. In all the RD cases, the energy is transferred from field 
to particle, and then radiation. Energy transfer, however, from field to particles becomes 
less efficient as radiation damping increases. Radiation damping suppresses the build-up of 
the high energy tail of the particle distribution. 

Figure El shows the momentum distribution of particles for (a) Run A(P X < 0) and Run 
C(P X > 0), and (b) Run D(P X < 0) and Run F(P X > 0). The v x B force accelerates elec- 
trons and positrons in the same direction along the x axis, whereas electric field accelerates 
electrons and positrons oppositely along the z direction, forming X shape distribution in 
the p x — p z plane as a result. The ponderomotive force J x B creates successive 'potential 
wells', which captures and accelerates co- moving particles in the x direction. We emphasize 
here that there is no charge separation in the x direction. Particle momenta in both x and 
z directions are radiated away in both weak and strong magnetic field RD cases. Radiation 
damping is more severe for higher B. We also find that radiation enhances the bifurcation 
instability discussed by Liang and Nishimura^J]. 

Next, we compare the radiation power of RD cases with NRD cases. For the NRD case, 
we estimate the radiation power using the relativistic dipole formula j3] 

where Fn and F± are the parallel and perpendicular components of the force with respect to 
the particle's velocity. The bracket () indicates that we take the average of all the particles 
located within 30Ad from the pulse front. We plot Runs A and D with k = 10~ 3 and 10~ 4 
to compare with the RD cases. 

For the RD (k > 0) cases, the radiation power is calculated using the formula 

(P) = \f rad -v\. (5) 

We compared the result of these two formulae for the RD cases, and they matched within 
the linewidth. Therefore, We only show the result of Eq. (JJJJ) here. 

Figure El shows the instantaneous radiation power (P) for all runs. The initial peak around 



tQ ce = 100 is due to synchro-cyclotron emission 



and the estimated radiation power for 



the NRD case quantitatively matches with the RD case. At later times (tQ ce > 1000), 



however, less power is radiated in the RD cases than NRD cases, because energetic particles 
are self-consistently decelerated by the radiation damping in the RD cases. 

Next, we calculate the self-consistent radiation field and its angular dependence directly. 
Intensity / and polarization IT of the radiation received by the observer located at x are 



given by 



and 
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J(n,r) = £Mret> (6) 
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where 

_ e n x [(w - &) x f3j] 
1 c (l-n-PfR ' [) 

is the radiated electric field from particle i located at r, to is a unit vector in the direction 
of x — r(r), (3 = v(r)/c, and /3 = d/3/dt. We assume that |aj| ^> |r| so that to is parallel 
to x. The summation in Eqs. © and (JJJ) is evaluated at the retarded time r — t — R/c, 
where R — \x — r\. We take r = when the pulse front reaches the observer. To specify 
the direction of the observer with respect to the x axis, we introduce 6 and <p as 

to = (cos#cos0, cos 6 1 sin 0, sin 6 1 ). (10) 

The time dependence of detected intensity logwl from each particle along = 9 = 
before the summation is shown as a contour plot in Fig. |U^a) and (c). Since the observer is 
in the positive x direction, radiation comes from the particles near the front in the positive 
x region. Since energetic particles are bouncing back and forth within the ponderomotive 
potential well and slower particles are dropped off from the well, the radiation duration is 
broadened. 

Sample intensity / is shown as solid lines in Fig. EJb) and (d), indicating the detected 
radiation peak width is At ~ 20O" 1 = 2\d/c for Run C and At ~ SOOfi" 1 = 8\d/c for Run 
D. This suggests that the radiation pulse width increases with B, but we emphasize that 
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actual detected pulse are much broader due to angular effect. Polarization II is shown as 
dotted lines in Fig. E[b) and (d), confirming that the radiation is strongly linear-polarized 
as anticipated from the initial uniform magnetic field with the small depolarization coming 
from the initial random velocity distribution in the y direction. 

In Fig. 03 we show the total angular fluence J J(0, 9)dr as a function of the view angle. 
Angular fluence peaks around # = 3~8°in0 = O cases (solid and dash-dotted lines), 
corresponding to not the direction of the Poynting vector but the direction of high energy 
particles in Fig. |21 Angular fluence rapidly decreases with in 9 = case, indicating 
radiation is strongly collimated in the x — z plane. 

In summary, we observe the self-consistent radiation damping and radiation output from 
the acceleration of electron-positron plasma via a relativistic PIC simulation. We find that 
the coupling between the field and particles becomes less efficient with larger radiation 
damping, implying the actual cooling time by the radiation increases. Radiation damping 
decelerates the energetic particles, but the ratio between P x and P z retains. The resulting 
radiation angular fluence peaks in the direction of the energetic particles. The radiation 
power monotonically decreases in the early time and becomes more or less constant later, 
and it is always smaller than the non-radiative case, indicating continuous cooling by the 
radiation damping. The radiation field is strongly linear-polarized both in weak and strong 
magnetic field cases, which may be detectable by 7-ray burst observations as an indication 
of Poynting flux acceleration. The simulations shown here are still too short to determine 
the cooling time when most of the EM energy is radiated away. Such questions remain to 
be answered by much longer runs. 

This research is partially supported by NASA Grant No. NAG5-9223, NSF Grant No. 
AST0406882, and LLNL contract nos. B528326 and B541027. KN is also grateful to ILSA, 
LLNL for hospitality, and B. Remington and S. Wilks for useful discussions. 
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TABLE I: The parameters for each runs 



k ujpe/Qce Duration tfl ce 



Run A 0.1 10000 

Run B 10" 4 0.1 10000 

Run C 10~ 3 0.1 10000 

RunD 0.01 70000 

RunE 10" 4 0.01 70000 

Run F 10" 3 0.01 70000 
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FIG. 1: System-integrated energy in the electromagnetic field (1), particles (2), radiation loss (3), 
sum of field and particle energy (4), and total energy (5) as functions of time for Run D (solid), E 
(dashed) and F (dash-dot). 



11 



£ 




-25 



-50 
-300 -150 



p /mc 



a 
E 



150 300 




-250 



-500 



-5000 -2500 



p /mc 



2500 5000 



FIG. 2: Comparison of momentum distribution of particles between Run A [(a), P x < 0] and Run 
C [(a), P x > 0] at tVL ce = 5000 (circle, asterisk) and tQ, ce = 10000 (x, plus), and between Run .D 
[(b), P x < 0] and Run F [(b), P x > 0] at tQ ce = 35000 (circle, asterisk) and tfl ce = 70000 (x, plus). 
Results for P x > and P x < are identical in all cases, so only half momentum spaces are plotted. 
This shows that radiation damps both P x and P z but the ratio between them relatively unaffected, 
and that the damping increases with magnetic field. 
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FIG. 3: Average instantaneous radiation power from particles within 30Xd from pulse front for 
(a) tu pe /{l ce = 0.1 (Runs B, C) and (b) 0.01 (Runs E, F). Estimated power for NR cases (Runs A, 
D) are also shown for comparison. Radiation monotonically decreases in time, and reaches more 
or less constant in (b), whereas it still decreases in (a). 
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FIG. 4: Contour plots of instantaneous intensity log 10 1(9 = (p = 0) as functions of source 
frame coordinates (x, t) for for (a) Run C and (c) Run F. Intensity (dashed lines, left scales) and 
polarization (solid lines, right scales) as functions of detector time r for (b) Run C and (d) Run F 
at 9 = <j) = 0. Intensity scale is arbitrary. To obtain time dependence of intensity, instantaneous 
intensity from each particle is summed up along the light cone r = t — R/c =const., shown as a 
black dashed line in panel (a). Note that actual detected pulses are much broader when we take 
into account angular effects (c.f. Fig. 03) 



14 



10 



20 



10" 



"i 1 1 1 1 1 r 



RunC, <|>=0 

RunC, 9=0° 

Run F, <j)=0° 




J I I I I I l L 



2 4 6 



8 10 12 
<|), 6 (degree) 



14 16 18 20 



FIG. 5: The total radiation angular fluence J I(cp,9)dT for Runs C and F as a function of the 
angle (</>, 9). 
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